Anisotropic low-temperature piezoresistance in (311)A GaAs two-dimensional holes 
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We report low-temperature resistance measurements in a modulation-doped, (311)A GaAs two- 
dimensional hole system as a function of applied in-plane strain. The data reveal a strong but 
anisotropic piezoresistance whose magnitude depends on the density as well as the direction along 
which the resistance is measured. At a density of 1.6 x lO'^^ cm~'^ and for a strain of about 2 x 10~* 
apphed along [Oil], e.g., the resistance measured along this direction changes by nearly a factor 
of two while the resistance change in the [233] direction is less than 10% and has the opposite 
sign. Our accurate energy band calculations indicate a pronounced and anisotropic deformation 
of the heavy-hole dispersion with strain, qualitatively consistent with the experimental data. The 
extremely anisotropic magnitude of the piezoresistance, however, lacks a quantitative explanation. 
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The piezoresistance effect ~ the change of a device's 
electrical resistance as a function of applied stress - is of 
great technological importance as it is utilized in making 
force, displacement and pressure sensors [jy]. It has also 
been used recently in measuring sub-nanometer displace- 
ments of tips in atomic force microscopes 12^, ^ , and pro- 
posed data storage devices based on piezoresistive read- 
back [J]. 

Strain affects the resistance of a solid in two ways. It 
changes the physical dimensions of the resistor and can 
also modify its energy band structure. In metals, it is the 
former effect that is dominant and the the gauge factor, 
defined as the fractional change in resistance divided by 
the solid's fractional change in length, is dictated by the 
Poisson's ratio which is typically around 2. On the other 
hand, for certain semiconductors the modification of the 
band structure is a more prominent factor affecting the 
piezoresistance. In multi-valley semiconductor systems, 
for example, strain changes the relative population of the 
valleys. If the valleys have anisotropic effective masses, 
the charge transfer induces a large resistance change in 
certain crystallographic directions Gauge factors of 
~1500 have been reported for n-type Si inversion layers 
and ~10,000 in AlAs two-dimensional (2D) electron 
systems at low temperatures |3] . 

In this Letter, we report a large, density (p) dependent 
piezoresistance effect in a (311)A GaAs 2D hole system 
(2DHS). Wc utilize a simple but powerful technique to 
apply quantitatively measurable in-plane strain in-situ 
Q and measure resistance parallel and perpendicular to 
the applied strain. We deduce a gauge factor of ~ 3600 
at p = 1.6 X 10" cm-2 and T = 0.3 K. The data can be 
qualitatively understood by the change in the heavy-hole 
valence band shape with applied strain. But a quantita- 
tive explanation of the observed piezoresistance magni- 
tude and particularly its extreme anisotropy, is lacking. 

Our sample is grown on a (311)A GaAs substrate 
by molecular beam epitaxy and contains a modulation- 
doped 2DHS confined to a GaAs/AlGaAs heterostruc- 
ture. The Alo.35Gao.65As/GaAs interface is separated 



from a 17 nm-thick Si-doped Alo.35Gao.65As layer (Si con- 
centration of 4 X 10^^ cm~'^) by a 30 nm Alo.35Gao.65As 
spacer layer. We fabricated L-shaped Hall bar samples 
via photo-lithography and used In:Zn alloyed at 440°C 
for the ohmic contacts. Metal gates were deposited 
on the sample's front (lOnm Ti; 30nm Au) and back 
(lOOnm Ti; 30nm Au) to control the 2D hole density. 
Magneto-resistance data taken at p = 2.1 x 10^^ cm^^ 
and T = 0.3 K yield mobilities of 1.7 x 10^ cm^/Vs and 
4.3 X 10^ cm^/Vs in the [Oil] and [233] directions respec- 
tively. This mobility anisotropy in the (311)A growth 
direction stems from the quasi-periodic ridges along the 
[233] direction [§,[l|. 
We apply tunable strain to the sample (thinned to 
200 ^m) by gluing it on one side of a commercial mul- 
tilayer piezoelectric (piezo) actuator with the sample's 
[Oil] crystal direction aligned with the poling direction 
of the piezo (Fig. [1] inset) [81]. When bias Vp is applied to 
the piezo, it expands (shrinks) along the [Oil] for Vp > 
{Vp < 0) and shrinks (expands) along the [233] direction. 
We have confirmed that this deformation is fully trans- 
mitted to the sample, and using metal strain gauges glued 
to the opposite side of the piezo, have measured its mag- 
nitude 8, llj. Based on our calibrations of similar piezo 
actuators, we estimate a strain of 3.8 x 10^*" V^^ along 
the poling direction. In the perpendicular direction, the 
strain is approximately —0.38 times the strain in the pol- 
ing direction [^IT3]. In this paper we specify strain values 
along the poling direction; we can achieve a strain range 
of about 2.3 x 10"^ by applying -300 < Vp < 300 V 
to the piezo Ts*]. Finally, the back-gate on the sample 
is kept at a constant voltage (0 V) throughout the mea- 
surements to shield the 2DHS from the electric field of 
the piezo actuator. 

Figure [1] summarizes our measured piezoresistance. 
Three trends in the data are clear. First, the change 
in resistance (i?) with strain is density dependent and 
is larger for lower densities. Second, the piezoresistance 
is more prominent along the [Oil] direction. Indeed, at 
p = 1.6 X 10^^ cm^'^, R along this direction changes by 
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FIG. 1: (Color on-line) Piezoresistance measured along the 
[Oil] and [233] directions for three densities. To obtain the 
resistivity (in units of kQ per square), multiply 7? by 0.3. 
Inset: Experimental setup. The poling direction for the piezo 
is along [Oil]. The strain values indicated in the top scale 
have an uncertainty in the absolute value and hence specify 
the change in applied strain (see Ref. [U)- 



a remarkable factor of two for an applied strain of only 
^ 2 X 10^^. Third, in the entire strain range we have 
studied, R along [Oil] is always larger than in the [233] 
direction. 

In order to understand the data of Fig.[T] we performed 
self-consistent calculations of the Fermi contours of the 
2D holes residing in the heavy-hole (HH) valence band 
as a function of strain. Note that in this system the spin 
degeneracy of the HH band is^ lifted due to the presence 
of strong spin-orbit coupling Jl4]. The spin subbands are 
termed heavy HH (HHh) and light HH (HHl), reflecting 
the magnitude of their effective masses. We used the 8x8 
Kane Hamiltoniari, augrnented by the strain Hamiltonian 
of Bir and Pikus flj, ll3 j to take into account the spin- 
orbit coupling of the system as well as the strain-induced 
contributions. We adapted this model to the (311) ori- 
entation of our sample by a suitable coordinate trans- 
formation. We note that these energy band calculations 
are quite accurate and have quantitatively explained the 
spin-orbit induced spin-splitting and its dependence on 
density, electric field and strain in 2DHSs [igI. [TtL [ll] . 

The results of the energy band calculations for our 
sample parameters are shown in Fig. [2] Since all the 
densities show similar trend, we only discuss results for 
p = 1.6 X 10^^ cm~^. Figure [2] shows that the anisotropy 
of the HHh band is strongly enhanced with the appli- 



FIG. 2: (Color on-line) Calculated Fermi contour plots for 
the HHl and HHh spin subbands; A:[oii] and fcpsaj are the 
wavevectors along the two perpendicular ([Oil] and [233]) di- 
rections. Inset: Comparison between the experimentally mea- 
sured (symbols) and calculated (lines) piezoresistance. 



cation of strain. This results in the anisotropy of the 
HHh effective mass as well. Qualitatively, a larger Fermi 
wavevector in a particular direction typically implies a 
larger mass in the same direction. Hence, applying ten- 
sile strain along the [Oil] direction leads to a larger mass 
along [Oil] and a smaller mass along [233]. This change 
in effective mass qualitatively explains the change in re- 
sistance in the two directions. 

To make a direct comparison with the experimental 
data, we calculated the resistance for our sample param- 
eters taking the strain dependence of the anisotropy into 
account. We assume a constant and isotropic relaxation 
time (t) and sum over spin subbands, n, with energy dis- 
persion En(k). We evaluate the velocity, dEnJ (hdk) at 
Ep and T = to obtain the conductance (a) pjj] , 



(Pk dEn dEn 



(27r)2 hdh hdk:, 



S[EF-En{k)]. (1) 



Since it is difhcult to estimate the value of r, we normal- 
ize the calculated resistance to its value at zero strain 
and compare the calculations to the data of Fig. [TJ For 
the comparison shown in the inset of Fig. [2l the exper- 
imental resistance was also normalized to the resistance 
at zero applied piezo bias, Rq. As the figure shows, the 
calculations certainly match the qualitative trend seen in 
the data. It should be noted that the calculations are 
based on the particular sample structure and density but 
contain no fitting parameters. 

A puzzling aspect of Fig. [5] (inset) data is that in 
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FIG. 3: (Color on-line) Electric field dependence of piezore- 
sistance in a 20 nm square well sample grown on a (311)A 
surface at constant density. The gate biases are: (Vf ~ - 
0.2 V, Vi3 = V) for Ej_ = -4.5 kV/cm, (Vf = -0.375 V, Vs 
= 84 V) for E± = 0, and (Vf = -0.5 V, Vs = 144 V) for Ei_ 
= 3.2 kV/cm. 



the [Oil] direction, the experimental change in resistance 
is typically larger than predicted by the calculations. 
The [233] direction, on the other hand, exhibits a much 
smaller change in the measured resistance than in the 
calculations. It is tempting to relate this mismatch to 
how the quasi-periodic corrugations [9j parallel to [233] 
might affect the directional dependence of t with strain. 
However, our piezoresistance measurements on a 20 nm 
square quantum well sample, as a function of perpen- 
dicular electric field {E±) at a constant density, indicate 
otherwise [Fig. [3]. The square well sample is glued to 
the piezo with the poling direction in the [Oil] crystal di- 
rection. E± is tuned while keeping the density constant 
with the help of front (Vp) and back (Vb) gate biases 
[20| . At E± = 0, the confining potential of the square 
well and the carrier wavefunction are symmetric. For 
E± 7^ 0, the potential becomes asymmetric, resulting in 
the carriers being 'pushed' closer to the interface. Hence 



the resistance increases with the increase in the magni- 
tude of E± because of interface roughness, as shown in 
Fig. [3] at a particular piezo bias. If the direction de- 
pendent mismatch of the piezoresistance were related to 
the corrugations along [233], we would expect a larger 
piezoresistance mismatch with an increase in the mag- 
nitude of E± . On the contrary. Fig. [3] shows that the 
piezoresistance is nearly independent of E± [21]. 

Moreover, a similar mismatch was also observed in 
(100) GaAs 2DHSs [l^l, where the change in resistance 
was reported to be larger than the calculated values 
in the [Oil] direction, while along [Oil] the measured 
change was much smaller than predicted by the calcula- 
tions. The mobility anisotropy in (100) GaAs 2D holes 
at zero strain is only about 1.2-1.3 between the two per- 
pendicular ([Oil] and [Oil]) in-plane directions. This is 
much smaller than the (311)A samples where the mo- 
bility anisotropy is typically about 2 to 3 between the 
[Oil] and [233] directions. It is unlikely, therefore, that 
the corrugations in (311) A samples are responsible for 
the mismatch between the rate of change of resistance 
along the [Oil] and [233] directions. We emphasize that 
this mismatch is independent of the poling direction of 
the applied strain: Our measurements on a different het- 
erostructure sample from the same wafer as the one used 
in Fig. [1] but glued on a piezo actuator with the poling 
direction along the [233] direction also exhibit a larger 
change in resistance in the [Oil] direction compared to 
[233]. 

In conclusion, we observe a large piezoresistance ef- 
fect in (311)A GaAs 2DHSs. The data indicate a strong 
dependence on density, and on the in-plane direction 
along which the resistance is measured. The maxi- 
mum strain gauge factor we measure (along [Oil] and 
for p = 1.6 X 10^^ cm~^) is ~ 3600. This value is com- 
parable to the ones reported, at liquid He temperatures, 
for p-type Si (110) inversion layers at p = 5 x 10^^ cm~^ 
[11] and for GaAs 2DHS grown in the (100) direction at 
p = 7.6 X 10" cm-2 22]. 

We thank the ARO, DOE and NSF for support, and 
M. Grayson for stimulating discussions. 



[1] A. Belu-Marian, E. Candet, and A. Devenyi, Thin film 
resistive sensors (TOP Publishing, 1992), p. 113. 

[2] M. Tortonese, R. C. Barrett, and C. F. Quate, Appl. 
Phys. Lett. 62, 834 (1993). 

[3] J. A. Harley and T. W. Kenny, Appl. Phys. Lett. 75, 289 
(1999). 

[4] B. W. Chui, T. D. Stowe, T. W. Kenny, H. J. Mamin, 
B. D. Terris, and D. Rugar, Appl. Phys. Lett. 69, 2767 
(1996). 

[5] C. S. Smith, Phys. Rev. 94, 42 (1954). 

[6] G. Dorda, 1. Eisele, and E. Preuss, Solid State Comm. 

11, 1625 (1972). 
[7] Y. P. Shkolnikov, K. Vakili, E. P. De Poortere, and 

M. Shayegan, Appl. Phys. Lett. 85, 3766 (2004). 



[8] M. Shayegan, K. Karrai, Y. P. Shkolnikov, K. Vakih, 
E. P. De Poortere, and S. Manus, Appl. Phys. Lett. 83, 
5235 (2003). 

[9] M. Wassermeier, J. Sudijono, M. D. Johnson, K. T. Le- 
ung, B. G. Orr, L. Daweritz, and K. Ploog, Phys. Rev. 
B 51, 14721 (1995). 

[10] J. J. Heremans, M. B. Santos, K. Hirakawa, and 
M. Shayegan, J. of Appl. Phys. 76, 1980 (1994). 

[11] O. Gunawan, Y. P. Shkolnikov, K. Vakih, T. Gokmen, 
E. P. De Poortere, and M. Shayegan, Phys. Rev. Lett. 
97, 186404 (2006). 

[12] The strain in the perpendicular direction is taken into 
account in our calculations described later in the paper. 

[13] In our experiments, we can determine the relative 



4 



changes in values of strain accurately, but not the piozo- 
bias corresponding to zero-strain precisely. This is be- 
cause, thanks to a mismatch between the thermal ex- 
pansion coefficients of GaAs and the piezo-stack, at low 
temperatures the sample can be under finite strain even 
at Vp = 0. This residual strain is cooldown-dependent 
and we do not know its precise value for our data. Based 
on our experience with cooldowns of samples glued to 
similar piezo-stacks, we expect a residual strain up to 
about ±1 X 10"''. 
[14] R. Winkler, Spin-Orbit Coupling Effects in Two- 
Dimensional Electron and Hole Systems (Springer, 
Berlin, 2003). 

[15] G. L. Bir and G. E. Pikus, Symmetry and strain-induced 
effects in semiconductors (Wiley, New York, 1974). 

[16] J. P. Lu, J. B. Yau, S. P. Shukla, M. Shayegan, 
L. Wissinger, U. Rossler, and R. Winkler, Phys. Rev. 
Lett. 81, 1282 (1998). 

[17] S. J. Papadakis, E. P. De Poortere, M. Shayegan, and 



R. Winkler, Physica E 9, 31 (2001). 

[18] B. Habib, ,J. Shabani, E. P. Dc Poortere, M. Shayegan, 
and R. Winkler, Phys. Rev. B 75, 153304 (2007). 

[19] R. Chambers, Electrons in Metals and Semiconductors 
(Chapman & Hall, London, 1990). 

[20] For details of the square quantum well sample structure 
and tuning the wavefunction symmetry see: S. J. Pa- 
padakis, E. P. De Poortere, M. Shayegan, and R. Win- 
kler, Physica E 9, 31 (2001). 

[21] This independence of the piezoresistance on E± is cor- 
roborated by our calculations of the Fermi contours as a 
function of E± for our quantum well. 

[22] K. L Kolokolov, A. M. Savin, S. D. Bencslavski, N. Y. 
Minina, and O. P. Hansen, Phys. Rev. B 59, 7537 (1999). 

[23] G. Dorda, I. Eisclo, and E. Prcuss, in Proceedings of the 
11th Int. Conf. on Physics of Semicondcutors (Academic, 
New York, 1972), pp. 1468-1473. 



